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PREFACE 


The work described herein was done at the Research and Develop- 
ment Laboratories, Hughes Aircraft Company, Culver City, California, 
during the industrial term of the Slectronics Ingineering curriculun. 
My control system project was suggested by the Laboratories, and was 
carried out with their fullest cooperation as a part of a large scale 
study of airborne digital control systems. 

I shall assume that the reader is conversant with current terni- 
nolosy in the comouter field (2). What few terms are not in general 


use I will explain. The words "order" and "instruction" are synonymous. 
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I. INTRODUCTION 

This thesis describes in detail a orogram and coding for a 
particular aircraft navigation display control system utilizing a 
digital comouter, in less detail the equinment used, and finally a 
number of suggestions for programming of control systems. 

Some modern control systems become very complex and require 
many individual computing or comoensatineg elements. It seems likely 
that & gain might be achieved by combining all these comouting and 
compensating functions in one elaborate general-purpose comouter, 
perhaps on a time sharing basis (1). Hither digital or analog tech- 
nigues could be used. 

A digital computer offers certain inherent advantages in this 
aoplication. Chief among these are, first, that changes in function 
of parts of the control system, such as adjusting compensation of a 
servo loop, are made simply &s changes of numbers or instructions 
within the computer. No change of physical equipment is needed, 
Second, and more basic, is the memory, or storage of information 
that a digital computer contains. Thus, the control system can be 
influenced by its oast verformance, or by outside stimulation, Таке, 
for example, an automatic factory. Then the output of the factory 
could be influenced by inventories and sales records which had been 
gathered and collated by itself, 

Such an elaborate control problem is presented in seme high-o2er- 


formance aircraft, and the comouter on which the doscrived program was 
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run was a digital conouter for airborne use (3). 

Programming and coding a problem for à digital comouter have 
received much attention (4), but the &ooroach has bcen one of re- 
ducing to a minimum the human labor, either by admittedly inefficient 
computations (5,0¢.115) or by making the machine do the work of coding 
(6). This approach leads to overall efficiency when e oroblem is done 
once or a few times. If, however, a oroblem is repeated ad infinitum, 
as in a control system, human olanning time starts to nay off as re- 
duced ohysical equipment due to the lessened storage caoacity needed 
for the more efficient set of instructions. This will be discussed 
later. 

In any machine that is very comolicated some means of diagnostic 
testing should be provided. ‘this will be elaborated.. 

Finally, a reasonably didactic descriotion of the program will be 
given. ‘this orogram was a first attemot, and it is on the basis of the 


mistakes made that the recommendations have been formulated. 


(2) 
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CHAPTER I, DESCRIPTION OF EQUIPHINT 


The equipment used, while not the subject of this thesis, did 
influence the programming and coding, and so should be described. 
Peculiarities of the comouter logic important to the orogram will be 
described in later chapters when applicable. 

The conputer used was a first model which had been used for some 
time for circuit and component reliability and logical checks. No in- 
put-outout equipment had been designed soecifically for it, although 
input-output equioment was available that had been used on a very 
similar computer. One function of the Control Box (Chap. II) was to 
adapt and interconnect the computer and the input-output equioment. 
The Brown Teleplotter was used intact excent for changes in its in- 
ternal servo system to resoond to dc signals of the desired amplitude. 

The comouter is of the gensral purpose, serial, binary digital 
tyoe (3). It has an arithmetic unit, a control unit, and a magnetic 
drum memory unit. The arithmetic unit „erforms the aritnmetic oper- 
ations of addition, subtraction, multiplcation and division. The 
unit consists of three one work circulating registers and a binary 
adder, The magnetic drum memory (3) provides one eight word circulat- 
ing register, the "a" register, and fourteen sixty four word storage 
bands. These bands are interoreted by the control unit as eight order 
bands containing instructions, and six number bands. One of the number 
bands and one of the order bands are tied together, so that an instruct- 
ion can be interpreted as a number and changed, written back into the 


menory, and then used as a new instruction. Such changes are made as 


(3) 
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the result of the comoutations in progress, The instructions avail- 
able to the coder are shown in Table II. 

The input-outout equioment converts analog information (de volt- 
aces) into an equivalent binary number at a sampling rate of once every 
fourteen magnetic drum revolutions, writes the most recently samoled 
binary number onto the drum, reads the most recently written output 
binary number (written by the comouter) and converts it to an equiva- 
lent analog quantity (dc voltage). It is multichannel, having nine 
jnouts and four outputs, with one calibrating period (7). 

the plotting board consists of a chart and oen so arranged that 
two input dc voltages are interoreted and plotted as x and y deflect- 
ions of the pen on the chart, A potentiometer error detector feeds a 
modulator, whose ac output is proportional to the error. An amplificr 
then drives a motor in such a direction as to tend to make the error т го. 
This is duplicated for the two axes. The speed of resoonse of the oen 


is much faster than the contemplated olotting speed. 


(4) 
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CHAPTER II, TH PROBLE AND ITS PROGRAMING. 


1. The Problem 

A navigation display was desired to demonstrate the ability of 
the Hughes Model I digital comouter to store much information for a 
long time, The operation of the device wes visualized as follows. 
First, a map of the desired track is placed on a plotting board con- 
taining a movable plotting ocn. ‘The ven is placed in turn over the 
departure point and each of the enroute way points, or "fixes", At 
each point a button is depressed which will enter the oosition into 
the memory of the comouter. After all points are entered, the pen is 
returned to the departure point and the computations started by throw- 
ing a switch. The pen then "flys" to each of the fixes in turn at a 
fixed ground sneed on strä&csht line courses. After getting to the 
last fix, or destination, tke computer would turn itself off after 
preparing the program for entering the next set of fixes. Provision 
was made for entering sixteen fix positions. A signal lisht would in- 
dicate when the comouter was ready to have fixes entered, or ready to 
TUM. 

The Control Box had the following operational controls: two hand- 
wheels for East-West and lorth-South positioning of the oen during the 
entry portion of operation, a oushbutton for entering fixes and start- 
ing fly operation when ready, a function switch with Inter and Fly 


positions, and a switch to turn on recording pen on plotting board to 


(5) 
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prevent smearing. Contained in the same box were various inter- 
connections between the computer, the input-outout equipment, and 
the olotting board, various disconnect switches for the power and 
timing pulses, and reading and writing amolifiers required by the 
input-outout equipment. A schematic of this box is included as 
Pigure 1 to indicate the functions, although the details are im- 
material to this thesis. 

It was desired to make the fly routine a feedback looo so that 
the present position of the oen could be sensed in the computations, 
and so that smoothing could be included by placing a low pass filter 
after the dc analog outout. Smoothing could have been done in the 
program, and orobably will be included at a later date by Hughes 
engineers. Smoothing in the program would insure, for examole, that 
comoutational errors of a random nature would be smoothed out of the 
outout, rather than being cumulative. Smoothing in the dc analog out- 
put will accomplish the same result, but requires additional equioment. 
Snoothing should not be confused with the essential filter required to 
smooth the ina resbes 

In the course of the fly computations, the present position is 
written into the memory by the inout equipment. The comouter then 
reads this position and compares it with the fix position to which the 
"plane" is now headed. By a trigonometric method involving the fixed 


soced towards the fix position, an incremental change in position for 


(6) 
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the comouting period is determined and added to the present position 
to give a future position. This future position is then fed to the 
outout equipment, through the low pass filter, and to the pen to oosit- 
ion it. This future position then becomes the present position for 
the next computation. The comoutations are repeated every 14 drum 
revolutions ог 0.1 seconds, i.e., the samoling rate is ten cps. “ote 
that if a computational error is made and the plane gets off course, 
it will still continue to fly to the fix on the nearest straight line 
course. When the plane arrives at the desired fix, the comouter will 
change its instructions so that the next fix position will be used. 
After getting to the sixteenth fix oosition the computer will stoo, 
after setting in the proper instructions to enter new fixes. 
2. The Equations 
The course and incremental oosition are computed by the following 
set of equations. 
Let e y be the ——— of the present position. 

Хо,Уо be the coordinates of the next fix position. 

dt be the sampling period, 0.1 second. 

dx, dy be the coordinate distances travelled in time dt. 

У be the soecified ground speed towards the next fix. 


e be the angle of direction of travel, measured from 
Lhe, = axis. 


Then we can say that 


Yo "Y = tan 8 = ду/дх. 
X-X 


(7) 
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dx = ὙΠ cos 8. 

dy. GX #en S. 

A138 ¿asa 

Ynf 1* Ya, 4 dy. These last two are the future 
position, and become x, y for the next period. 

Difficulties arise from two sources. (1) Ambiguities due to 
sign of tan €, and (2), the fact that the computer can handle only 
numbers whose magnitude is less than one. These will be taken care of 
by using the conditional transfer ckS instruction. (Table Ша и ote 
that no instruction is provided to soecifically handle square roots, 
so an approximation will have to be made for cos 8. 

The equations finally used that follow, resolve these difficulties. 
See fig. 2. 

Given X, y, Хо» Vet, 2. 


1. Ріпа |5, -X]. 


2. Find |Yg -Y|. 
E wind 12,1 + (У рої. Store for future use. 
4, Find |Xy-X|l - |Yp-Y|. Check sign. I£ 
5a. Positive. Find 5b. Negative. find 
(Үр-Ү)/ 0-Х) = фәл Ө. (Xy-X)/(Yo-Y) : cot e. 
ба. Approximate cos 8 by polynomial. 6b. Approximate cos 8 by poly- 
nomial. 
Ta. Multioly Vàtcos9sz dx. 7b. Multiply Vdtcose9 = dx. 


(8) 
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Ва. Multiply dxtane@ = dy. 8b. Divide dx/coté = dy. 

2. Divide dx/(X,-X). Check sim. If 

10a. Positive. Do nothing. lOb. Negative. Subtract 
= de. 0 - dx = der, 

ll. Divide dy/(Yg-Y). Check sign. If 

l2a. Positive. Do nothing l2b. Negative. Subtract 


qa ut. Q.- dy =. дуй 


ዞ 


13. Add X / dx* z XQ,,... Write as output. 
14. Add Y # dy* < Yney =. Write as output. 

15. Subtract (|Mp-XI4(Yo-Y1) - z. Check sign. If 

lóa. Positive. Start over l6b. Negative. Start over 

with same Ap, Ур. with new Xp+1, Yps+ 1. 
17. Count wo a number such that after executing 16b sixteen times, 
a change of sign occurs to stop comouter. 

In the above, step 4 determines whether tan © is going to be 
greater than one. If so, cot 8 is used. Steps 9and 11 insure that 
dx and dy are of the proper sign. Steo 15 checks how closely the 
plane is to the fix position. hen this distance, or rather a dis- 
tance which is slightly greater than the distance to the fix, becomes 
less than some small number z, the olane has arrived, and the next fix 
position is used. 

The polynomial used to avoroximate cos 6 from tan 9 or cot 9 is 


derived by a Chebyshev volynomial operation, and is accurate to 0.3/9 


(9) 
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within tne range used: 


cose = 0.147 40. 4340x%40.9965, where х = tadê, or 


соз =у(.1474у ®-.4340у210.9965), where y = cote. 

Tnis approximation is easily done on the computer, and is easier than 
conventional algebraic methods would be. 
3. The Fly Program 

Using the instructions contained in Table II and these equations 
for & guide, the program was designed and is presented in Table IV. 
Table III shows a flow diagram of the instruction sequence. Before 
detailing the instructions some general explanations are required. 

On а machine which has available only such unsophisticated in- 
structions as shown in Table II, the detailed coding must proceed 
closely with the program,since the order in which operations are ver- 
formed depends upon, for example, when certain numbers become avail- 
able from the memory. It may be advantageous to delay some operations 
till a more convenient time. Thus the sequence of instructions may 
seem odd until the coding is studied in Chanter III. 

The subscript on instructions involving the D register denotes 
the sector number in this eight word memory. 

The B register acts as an accumulator for the arithmetic unit; 
that is, any number written into the B register adds algebraically 
to the number already there. In starting a new operation it is neces- 


sary to insure that the B register is cleared, as in instruction # 1. 


(10) 




















> æ š! Ú об ТВ нн ATI ገቴ 





что 4 “у 





ооа o» ww 2 b. о lel eo =e tim з к . 





е = ° tym» Wd 

= ee 

own? u в © σε Е m. MA 
Ij по 61 μα σε πμ... | 





> о. (ለመ ве === 'ኑ!ህ x. 
ÅJ a- Å I‘ j AN ame 
СЕЕ СЕТЕ ЕЕ ex (i 


— sal fw. se . ‹ ጋ 


dm mm m are pom ee m о ርረድ ተዯሇ ሙ 4 
rm ав Ай, um  — 

w ¿wy у= k. — >>" > | с э@ 
WE w 1: ο τι. 


Tr @ om di @ eee At = ! πω 


«8 8 ባህ κ a es eme © በበበክ ..፡፡ ... 









nn e 





3885 ‹<.፡ ብመ ιτ ο σα ο ο σι 
ДА Ал» - т — a | <> ' к= 4 "ων є cm 
c. ο τσι. a 111. | (5 4. 6 


ο. - α. AA A ITI “£ ΙΙ. —— | 





Any instruction involving "from B", such as ba, bm, etc., will clear 
the B register, 

А11 other registers and the memory will hold a number until a 
new number is written over it, renlacing the old number with the new 
number. 

Instruction #1 will be called simply 11, etc. 

Fl clears the B register. 

#2 transfers the fix Xy to the D register, This instruction will 
have to be changed as we progress from fix to fix. # 3 does same for 


Т, 


p° but will;he a seriesph of instmctiens,,asmm@llebe geen imeChapter Ill. 


#7 and #8 bring in the values of oresent position, x, 
and yu, from the memory. These numbers are then divided by a scale 
factor (7) to become the corrected x and y, and are stored in D (up 
to #17). 

#18 to 738 are manipulations to obtain K'4L', which is a rough 
measure of the distance to со to the fix, see Figure 2. Some inter- 
mediate results are stored in D for future use. 

The quantity K! L! can be a number from O to #2. This computer, 
which is designed to handle numbers between -1 and #1 (modulo 2), in- 
terprets a number between 1 and 2 as a negative number. Thus #39 cheds 
this fact, and if negative, indicating that K'/L' is greater than 1, 


* е е е е? a й, 
writes in its place an arbitrary positive number, Z. 


i4 Бо #49 form K'-L' as a measure of the octant of O, If in an 
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octant such that tan 9 is less than 1, the computations use tan 9. 
If tan © is greater than l, the computations use cot $. This avoids 
overflow on the following division. #50 senses this fact. 

#51 to #54 form tan 6, OR #73 to #76 form cot 9. 

Whichever form is used, #55 to #72 form either cos 90 or cos 6/cot 
9 by the approximation shown on page 10. 

#78 and #79 bring K'-L'! into the A register again to check sign. 
This was an easier operation to code than to store the actual decision, 
#50, previously made, i.e., the decision was made twice. 

#31 to #93 form dx and dy for the case when cot @ is used. The 
Sign of dx and dy may be wrong. 

#94 to #102 form ду and dx for the case when tan @ is used. Again, 
the sign may be wrong. 

#103 to #118 check that dx is in the same direction (has same sign) 
ES X coordinate distance to the fix. If wrong, the sign is changed by 
subtracting from zero. 

#119 to #137 does same for dy. Note again that a definite possi- 
bility for overflow exists in the two divisions, which are used to check 
parity of sign, when the present position aporoaches a fix. dx and dy 
must be larger than the actual change in present position due to the 
time lag effect of the smoothing filter. This mistake was corrected 
by utilizing a multiplication to check parity of sign in the actual 


machine, but was left in this program to illustrate that not only must 
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you think in the language of the computer, but also in terms of the 
overall machine considered as a control system, i.e., its analog as- 
pects. 

#138 to #149 form the future position and write this position 
into the memory in a position where the outout conversion device can 
read it. 

#150 to #154 comoare K'4ZL', or roughly distance to fix, with a 
small arbitrary number, z. When the oen has approached within z of 
the fix position #155 senses this change of sign. If not within z of 
fix, several waits are necessary so that recycling back to the #1 
start of the program will be after the 14th revolution of the magnetic 
drum. This is to insure synchronism with the sampling rate of the in- 
put-output equipment. 

If within z of fix, #160 to #165 take #2, add a number to the in- 
struction such that the "words to wait" (see Chaoter III) are changed 
to give a new Xp and Yp, and write the new #2 back in position so that 
it can be used the next time around, ach time a fix is changed a 
nunber is added to the count number such that after the sixteenth fix 
is passed the count number becomes greater than 1, overflows, and 
appears as a negative number to #175. 1179 then writes O in output ^, 
which is connected through the outout device to a signal lamp, puttim 


out the lamp. Later we will see that the lamo had been out on after the 
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lest fix was entered. #180 to #132 write in #2 of Table V, the 
entry program (see next section), so that the instructions will 
be correct for its initiation. 

#183 uses the same order as some number not zero so that 7134 
will stop the computer after getting to the sixteenth fix, or the 
destination. 

4, The Entry Program 

In Table V are the instructions used to enter fixes into the 
memory. 

T2 to 76 change the "words to wait" in the same manner es in 
previous section. Since this operation is done before using the 
order, the first "old am instruction" will be some number, meaning- 
less as an instruction, which when added to will be the desired "new 
am instruction". This will be discussed further in Chapter III. 

#7 to #16 alter oresent position, which in this case is the 
desired fix to be entered, by the same scale change as in previous 
section. 

#17 is an examole of orogramming influenced by the coding re- 
quirements. It vas necessary to enter this nunber into the same 
position in the memory as the "old am instruction" had been. This 
required sixty four word tines, and so the intervening time might 


as well be used for some more useful puroose. 


(14) 





a. i в 11» 9 Bu 
се “ ሚን i= = | 
я > ΤΙΝ. ΠΣ aa 
መ. IO nmn я: П 
>> > - 4 4 
т! я — e! 81 ' 
[^ we Y.” | a 
00 UC L EJ Е =: 
pb? = == ከ i ий» 
το. ° በጠ: 
ве ew à "NUS C А‘ 
. РУ, e ο = ‹ እይ а 
=» — зі яю» ማሜ чт m 
— 8) 3 во 2 
የሽ ብ. = b шъ ша ιν 
! n '1 ыа E B | 
E st (ወ 


ей, 





#21 makes use of the number written by #17 as an instruction, 
&nd enters the fix position into the memory. 

Since the instruction following #22 must occur immediately after 
writing the fix position (See Chapter III), and the fix positions are 
in separate places in the memory, 723 becomes à series of identical 
instructions, but with different "words to wait". This provides a 
means of determining when the sixteenth and last fix has been entered 
without utilizing à count number &s was done in the previous section. 
The lest of this series is called #28. 

For the first fifteen points, then, #23 to #27 puts an arbitrary 
positive number in to B, a zero into outout 4 for the signal light, 
and stoos the comouter ready for the next fix to be entered. 

#28 to #35 writes the instruction which becomes #2 of the fly 
program in its proper olace, writes a number which is nearly full scale 
into outout 4 to light signal light, writes a zero in the count number 
used in the fly brogram, and stops the computer ready for flying. 

5. The Test Progran 

The test program shown in Table I was included for diagnostic 
testing of the comouter components. It was of a tyve used previously 
for testing the Moel I computer, anà performs the aritmetic functions, 
various transfers, and utilizes all the number and order bands of the 


memory. It does these operations on all binary numbers fron 2 Uo liege 
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This program is not 





sxolained in detail since it is not the writers 
worx, and is immaterial to this control system, except in its ooer- 


ational use. 
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CHAPTER III. CCDI*G T€ PEOBLZSM 


Coding the program develo»ed in Chapter 11 requires a detailed 
knowledge of the computer and how instructions follow one another. 
Table I, contained in the envelope in the bac’: of the thesis, lists 
the coded problem. Programming and coding were done nearly simultane- 
ously. 

Three separate problems are coded in Table I: the fly program 
whose instructions are labelled Y, the entry program labelled £, and 
ше test program labelled T. Ene instructions in mach of the, three 
are serially numbered whenever oossible. Table III shows a flow dia- 
gram of instructions for the fly program. Reference should be made to 
fables IV and V for the sense behind the coding. 

The Model I comouter uses a rolotive address code (3) in which 
the address cf the next instruction is given by specifying the memory 
band in which the instruction occurs, and the number of word times till 
the new instruction is ready to be read. ‘The system is relative because 
the sector address of an instruction (or number) is never given, exceot 
for the first of a series of instructions, which always is in sector l 
in the selected menory band. (The p +A sector is the mtn word after the 
arbitrarily designated sector 64 in memory). 

An examole vill follow to illustrate this point. First, the "book- 


keeping" shorthand should be exolained. The instruction itself is follow- 
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ed by à series of three numbers in Table I. The first is the words 

to wait after reading the order, and before executins the order. The 
s number soecifies the number band of the memory involved if mean- 
ingful; thus, the instruction f does not irvolve the memory, but am 
does. The third number soecifies the nemory bend of the next instruci- 
ion. For examole, in sector 06, band 1, the instruction md, 2,4,1 

3T occurs. This means that the number in number band 4 should be trans- 
ferred to the D register, sector 1, after a 2 word wait, and that the 
next instruction, HT, will be found immediately after execution in 
band 1. The instruction was then read out during sector (word) 706, 
the computer did nothing during sectors #07 and 108, transferred the 
number found in number band 4, sector #09, to the D register during 
sector 709, and read instruction wT during sector #10. Note that the 
D register is eight words long, and arbitrarily sector 1 of the D 
register coincides with sector l cz the sixty four word memory. ‘Thus 
the D register repeats itself starting vith sectors 9, 17, etc. or the 
memory. 

The permissible words to wait before execution contained in any 
instruction can be anything from 0 to 16 words, except the wait in- 
struction (wt) which can have 0 to 32 words wait. The wait instruct- 
ion is executed immediately, or rather it would be ከ208 to say that 


execution corresponds to reading the next instruction. For examole, 
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check instruction $159? in sector 52, band 2, which feeds into in- 
struction FLF. 

Wote that a particular operation to be performed by an instruct- 
ion may be wrong. As long as the address of the next instructicn is 
correct tne computer will continue to cycle instructions blindly. 

This is of particular importance in trouble shooting. The usual ого- 
cedure is to male the instructions cycle, then worry about what mathe- 
matical results are obtained. 

The three programs all start in sector #1, but in different 
memory bands. The starting position is selected either by a set of 
buttons on the computer, or by the last instruction read by the con-- 
outer. Thus, vhen instruction fl149* is read and the computer stos 
ready to start on the entry orogram, the address of the next instruct- 
ion is band 3. When the computer is sterted again by oushbutton, the 
instruction #1E will be read first. 

It is by means of the variable words to wait portion of an in- 
struction that sixteen different fix positions are entered and taken 


out of the memory. However, in the actual number representing a particu- 


2 


lar instruction, the comoliment of the words to wait appears. Thus, 
the four digits reoresenting the words to wait will aopear oddly when 
added to in the manner shown. When all four digits are zero, this 


means 16 words to wait; adding one least significant digit means 8 
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words to wait, etc. See number band 2. 

Rather than taking up the detailed coding step by step, the 
interested reader may follow through as much of Table I as necessary 
to get the systom in mind, Compare Tables T, III aná IV with emphasis 
on the operation of the conditional transfer instruction ckÉ ana the 
variable words to wait, 

It must be obvious to anyone who has followed the coding in any 
detail that there is much wasted time and space here. The oroblem has 
become not one of electronics or control system design, but rather one 
for an ingenious bookkeeper and for numerical analysis. Time did not 
permit comoletion, If this control system were tc be used in an im- 
portant application, however, such effort would be essential to a good 


design. 
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CHAPTER IV. GHYERAL DISCUSSION AID RECOMMENDATIONS 


The test program used in this device was a good one in that ell 
functions of the comouter were tested. However, diagnostic indication 
of troubles encountered was very limited. A test should give a positive 
indication to a serviceman of the location of troubles. The instructions 
available in even as unsophisticated a comouter as the Model I should 
allow this by making more liberal use of the ckO instruction throughout 
the test orogram shown in Table I. Since a digital comouter is so lit- 
eral and human coding so liable to error, one must know whether the ma- 
chine is at fault or not before spending the hours required to trouble 
shoot a program for a control system. 

Once assured of proper machine operation, several steps can be 
taken to get a program to cycle instructions correctly. As noted pre- 
viously, it is usual to cycle instructions, then run test problems to 
check mathematical operations. After entering a set of instructions 
into the comouter in which a mistake is made, imoroper cycling will usu- 
ally result in circulation through an improper sequence of instructions. 
If this occurs, the following steos are of value with the Model I comouter. 

1. Change an instruction, say about half way through the orogran, 
to a ckO at some point where the computer will surely stop. If this 
half of the instructions cycles correctly and gives oroper intermediate 
results from a test problem, out the instruction back and change another 


instruction say three quarters of theway through. Continue this pro- 
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cedure till the trouble is narrowed down. 

2. After narrowing down, change one or more instructions to 
cycle only à small portion of the orogram reoeatedly, and get it to 
work. This might be called reentrant checking. 

3. Another method after narrowing down is to use test equioment 
built into computer (provided in most comouters) which will allow exe- 
cution of a single instruction at a tine. Check each sten in a suspect- 
ed region. 

Ц. In troublesome cases it may be necessary to check each in- 
struction with an oscilloscope. Аз might be suspected, this is time 
consuming and should be avoided if possible. 

5. A good deal of ingenuity helps more than set rules. 

When proper cycling is accorrolished some previously prepared test 
problems should be used to check operations. The method used with this 
set of programs was to disconnect the input equipment, write in known 
numbers where the inout numbers would have been written, and check either 
the final output, or, the numbers held in the registers at some intermedi- 
ate time, using the narrowing down technique previously described. 

If now the comoutations are orooerly made by the computer, the rest 
of the control system comes under scrutiny, especially the stability of 
the system and errors. Note that the program in use is an integral part 
of the control system, and must be thought of as part of the equioment. 


A useful concept in present control systems is the transfer function 
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either of the system or its components. A program may be thought of 

in terms of a transfer function involvins principally functions of the 
term exp(-as), a time lag (9). Unfortunately, too little has been done 
to develop a theory for this concept. It is mostly e study of numerical 
analysis (9), 

Damping of a control system has an analog in a digital comouter 
program. Derivatives and integrals are replaced in these discrete cal- 
culations with differences and sums. Damping then becomes smoothing, 
which may take the form of averaging or interpolation, or aooroximation 
of the required derivatives by differences. The problem is one of mmer- 
ical analysis and will be discussed further only in connection with noise. 

A simole method of estimating stability was suggested by Dr. Jacobi 
for the control system described. Since it is a feedback system involv- 
ing positive feedback, oscillation could occur when the total phase shift 
became a multiole of 2pi. The fly orogram has a time delay of 0.1 second, 
or a phase shift of f]= W/10, where м is the angular frequency. The RC 
filter used on the output equipment has a phase shift фо з WRC, where RC: 
0.5. Oscillation at w = O is meaningless. f = $4 j $2 z 2pi occurs 
when w = 15 pi, or at a frequency of 7.5 cos. However, the RC filter 
has an attenuation of 27 db.at this frequency so no trouble was expected 
or encountered. When the filter was disconnected for trial, oscillation, 
or at least random action did occur although the oen would of course not 


follow the gyrations of the numbers in the computer. Note that the same 
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smoothing effect could have been attained within the orogram had it 
been designed into the program. 

Errors within the comoutations (excluding a steady bias error) 
may arise from two principal sources. Round-off errors appear as 
white noise. large random errors may occur due to electrical inter- 
ference or errors in programming (see pg. 12), or due to marginally 
satisfactory equipment within the comouter. The random errors may be 
reduced by proper design and maintenance, but probably cannot be elinin- 
ated. Both types of noise can be reduced within the orogram by smooth- 
ing. White noise is most effectively reduced by some means of &verag- 
ing or interpolation. Large random errors can be treated, like ignit- 
Ton noise, by limiting the value of successive Increments in an iter- 
ative process, e.g., limiting dx and dy in the fly problem. 

An error which occurs as a bias error, €.£., a wrons scale factor, 
can cause considerable trouble in a system involving feedback. Scale 
factor is contained in the inout and output egvipments, within the pro- 
eran, and in any unas components external to these, particulariy the 
sensors. Thus the accuracy of the overall system is still determined 
by the weakest link. To take advantage of the possible accuracy of 
most digital comouters, then, attention should be directed to the de- 
sien of the analog components of the system, as well as to the program. 


In à time-sharing system it would be of value to utilize one computing 
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cycle for recomputing and writing a new scale factor into the orogram 
periodically, i.e., automatic calibration. 

Using a digital machine imolies both quantizing and sampling; these 
snould be consistent with system design. When hizh sampling rates are 
required, e.g., high roll rates of oresent jet aircraft, a stringent limi- 
tation is placed on the огосгат designer and the computer engineer to make 
an efficient (few words per cycle) and ravid (many words per second) ma- 
chine. Low quantizing error may be required in some aoplications. These 
two conditions are unfortunately contradictory in a serial machine. This 
is especially unfortunate when one equipment is used on a time-sharing 
basis for many jobs. The use of a parallel machine may be required with 
its higher speed of operation but greater amount of equipment, if the 
program designer is not dexterous enough. 

A program may be made more efficient by efforts of the numerical 
analyst in reducing the steos required to do the series of operations. 
The coding of a particular orogram is a related but different oroblen. 

In coding, each weit for the aopearance of a number in the memory із іп- 

efficient, as are comoletely filled sectors when other sectors are rela- 

tively vacant (compare sectors 5 aná 55 of Table I). The requirements of 
coding will orobably determine the details of the computer design in the 

long run. Efficient coding is the result of attention paid to a great 


many small details. 
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Figure lI. GEOMETRY OF FLY PROBLEM. 
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- 
Words to 
Symbol Meaning Zxecnte 
en Jpenster number in c register to specified 
| memory band l 
am Transfer number in a register to specified 
number band 1 
dm. Transfer number in d register to specified 
nunber band 1 
bm Transfer number in b register to soecified 
mzrnber band. b cleared. 1 
ma Transfer number in specified number band to 
a register 1 
mc Transfer number in soecified number band to - 
c register 14 
ша Transfer number in soecified number band to 
d register 1 
ca, Transfer number in c register to a 1 
Aca Transfer number, disregarding sign, from с 
register to a I L 
da Transfer number in d register to a 1 
ba Transfer number in b register to a register. 
b cleared 1 
σα Transfer number in c register to d. 1. 
ад Transfer number in a register to d. 1 
de Transfer number in d register to c. 1 
ра Transfer number in b register to d. l 
b register cleared 
+ Add number in a to number in b, sum in D 
register 1 
- Subtract numser in a from number in b, 
difference in b 1 
X Multiply a by c, result in » resister 16 
div Divide Ὁ by, resukts in-cyeb=aldeaned <q 16 
wt Do nothing for (1 to 32) words 0 
ско Check b register. If zero go to next order. 
En Ut, stos9*"osubuter 3. 
cks Check sign of number in a. 1f positive - - 1 
If negative - - 2 


Table II. Instractions avuilable in Model I 
Digital Comouter. 
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Problem 


sign of IXp-X1 2 |Yy-Yl 


sign of |Хо-Х| =. 1፡ --ቼ] 


direction of dx 


direction of dy 


arrival at fix 
up words to wait of order 
Sepial number of fix 


159 375, 
176” “77 Prepare for entry routine 


184 Check zero. Computer stops 


Table III. Flow Diagram of Fly Instructions 


28 


meldouis 


Y-a% X X-gX to mela 


5-95 š Χ-αΧ Ίο በኋፈ8 


Y-gY - X-gX “To በ228 


хр To noljoerib 


vb 1ο поідовчіб 


х13 Je вт: ччв 
чорто To Фіви 02 аБтон qu 
Kit to tedmun [sites 


18436 


Xoeno — 


As 8በ5 


Hoodoo 


Ho edo 


Ho endo 


Хоейо 
Зло 
Хоейо 


еп г иоч vine чо? ет8ае9ч4 УУД 


в4о?в чефицшо0 „отед Хоейо ASL 


впо| дома вит ҮЧ Зо швчдв! 4 wol" .III eldst 


8S 


ፈቱ 


єв 


oy 


£8 
ce 


oL I 
SLL 


C 


PEL 


004 
СТІ 


ec 


EA 


га 
E. 
аа 
4ο 
SOL 
EOL 
ALL 
aL 


ell 


ዕርፊ 
ІСІ 


Sti 
Qu 


сет 
ací 


осі 
ат. 


ва M m m r s тИ п r 


Bo am C Reg. e. in "йогу 





0 
Ὃν 
їр 
Yo in Da 
А in Ds 
Ay in 24 
τα 
Yu 
6. 
0 . 
Y dm D 
Ku 
Ха 
S.C. 
0 Χ 
in De 
fp 
Mb 
X 
Хот 5 5 
0 t in D, 
K 
"як 
E" din D 


Table IV. Program of Fly Instructions. 
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19 „ал K'=L! in m. 
O. ONS 8. ጩ to 51 
Fee 0 73 
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го гре. г Пи 








insts 


51.Аса 
52.4 
22. 925 
54.div 
55,65 
56,591 
57.Х 
58. be. 
59.ado 
60,92፡ 
61.X 
62.55 
63.ma 
6!ν,Χ 
65,bd4 
66.dsc 
67 „ша 
68,X 
69. ἄπο, 
70.4 
7i.ma 
Vlas 
73.495. 
74.4 
79 Sca 


lo.in A der. No.in B her, To.in C Per, 





+ 
IM 
ті 
K 
0 L! /K<=+*tanə 
tan9,cot9 
tan“ 
tan“ 0 
tan? 
tan! 
0 
E 
-. 43 бад“ 
0 
tan? 
0.1474 
„147 ел 
-„ЦЗЦфол? 
sum 
0.9965 
next inst 77 cos 9 
Kt 
Kt 
Τι; 
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о. YA menor 


Jj. 7 


tan in Dy 


tan^ in D^ 
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tam in Ds 


in Dy 


$ 
б 

І, 
—_ с መሙ 

е 

"a. 
a. 

= су 
"1. 

c8 t 
— . 

.. 
"πι 
«ፍ: =j a 


в аз «І slim 





Мо. № А Ней. doin В Res, Хо, 12 С,Веє, Jo. in memory 


. 
76.div Next inst 55 0 K' /L!z coto 








77.00) 0 cos or cos/cot 
in D4 
79.wt 
79.ma Kt-L! 
80.ckS if +, to 94 
if το Dl 
61,966 cos/cot 
02.419 cot 
II cos 9 
85, bdj, 0 cos in Dh 
85.шс Vdt 
86.dıa cos 8 
87.2 a 
88.ba dx 0 
89,502 dx in До 
90.4 dx 
91.d1& cot © 
92.div 0 dy 
93.сд5 next ins 103 dy in ds. 
9H. uc Vdt 
95.аца, cos @ 
96.X dx 
97 „Ба, dx 0 
93.ad5 dx in Do 
99 .41с tan 9 


100.X dy 
Table IV. Continued 
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4 ite of) £u 


= sp 


ет | ро =. mel 


No.in B lI Nosin C Row. o, in memory 





-101.bds 0 dy in Dg 
102.wt 

103.42 X. 

104.4 Хр 


105. дса x 





0. Xy AE 
107 .bdy, 0 Kin Dy 
108. doa dx 

109. + ах 

110.ада K 

111.div 0 dx/X 

112.cd: ax/X in D] 


113.496 dx/K 


11^,ckS i$ £, e 115 
if ~, to 116 
115.wt next ins 119 
116. ἄρθ ጩ 
117.- -üx- dx* 
118.042 0 dx* in D, 
119.dga Ya 
120. Yo 


123., 925 Y 

122.- Yp-Ys1 

123.bdj 0 L in D4 
124. dca dy 

125.7 dy 
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=ч ~= 








Insi. 


126; аца 
lav 
128. са 
129.048 
130. ck5 
334.5 
132.7 
122.658 
135.- 
135. ba 
136.ад с 
137.7 
135, 525 
139.# 
140.545 
141. .д2а, 
142.4 
143. ρα 
144.4 
145,02. 
466 , 9.55 
147 .am 
148.241 
149. сп 


].50 „ва 


io.in À Rees. 


dy /L 


i$ й τῷ 131 
17 -, 50 133 
ду 


next ins135 


dy 


dx* 


x 


КИЙ! ог £ 


~Oein B Бег. lo.in C nes, 


0 ey /L 


dy 


-ду=ду* 


dy* 


ту дуката. 
0 


Жена. 
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Loin lismory 


dy/L in Ty 


dy* in Ὡς 


Y,71 in D5 


Yıurl in output 


Xnfl in Di 


Та 1 in outout 


AA - 





inst. kosin A Reg, No.in 3 Reg. ko.in Q Reġ. lip,in memory 





. 
151.4 Kt/Lt or £ 
152 .ma, 2 
159.- (X! ZL!)-z 
154. ba (፪፤#:,፤)=2 0 
196.055 if # to 156 
if -, to 160 
156.wt 
157.wt 
158.wt 
159.wt next inst 1 
160.wt 
161.ma old ma 
insimıctien 
162.4 old ma 
instruction 
163.ma 1/16 
164.4 new ma 
instruction 
165.bdy, 0 new ma inst 
Ada за 
166 „та, a 
መ) = 
167.4 27742715 
168.ma count no. 
169.4 new count no. 
170.ba new count no 0 
l7lwwt 
172.d4c new na 
instruction 
175. ci new ma inst 
in memory. 
174, ል: new count по. 


in menory. 
175.685 if Z,to 176 
зво I 
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Inst. Doin A keg, Mo.in B eg. Wo,in С Бек. llo, im memory 


- 
176.wt next insl59 
(recycles) 
177 .ma 0 
178.wt 
179.am O in output À 


130.ma new am inst 


131 .wt 

192.am new am inst 
in memory 

183.7 new am inst 


1 
18k, c10 stop comouter 
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My ^! 





Inst. Mo.in A Ros, 110 ,ኋቬ В Ree, Po im G Regs lo.in memorv 


* 
1.ра 0 
2. па old am 
instruetion 
3.4 old am 
anerruction 
4.ma 1/16 
5.7 new am 
і instruction 
6.54, 0 new am inst 
in Dy 

7 . ma, Xu 

3.mdg Y, За Dg 

9.4 x6 
10.ma ደ 
IN div 0 X 
12.452, Yu 
13.с4> X in Do 
14.4 Е. 
1 22486 ዑ. 
ТОБ 0 Y 

18. dm new am inst 

in memory 

19.wt 
20 .wt 


Table V, Program of Entry Instructions. 
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ET I 4 









105%. Lo,in ес, loin D Reg. losin 6 nes, o. in memory 


21.2m X in memory 
22 eî Y in memory 
λας 0 

21 „ша, 5 

25.7 3 

26.cm O in output + 


27.ck0 Computer 
stops 
23.me ma instr 


29.ma ie’ 


30.4 iuge? 

31 дещ ma inst in 
memory 

32.110 0 

33.am full scale in 
output 4 

35,6ዉ O in count no. 


35.ck0 Computer 
stops 


Table V, Continued 
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